Microscopic expression for the heat in the adiabatic basis. 
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We derive a microscopic expression for the instantaneous diagonal elements of the density matrix 
Pnn{t) in the adiabatic basis for an arbitrary time dependent process in a closed Hamiltonian system. 
If the initial density matrix is stationary (diagonal) then this expression contains only squares of 
absolute values of matrix elements of the evolution operator, which can be interpreted as transition 
probabilities. We then derive the microscopic expression for the heat defined as the energy generated 
due to transitions between instantaneous energy levels. If the initial density matrix is passive (diag- 
onal with p nn (0) monotonically decreasing with energy) then the heat is non-negative in agreement 
with basic expectations of thermodynamics. Our findings also can be used for systematic expansion 
of various observables around the adiabatic limit. 



The standard definition of heat is the energy trans- 
ferred from one body to another due to a temperature 
difference. However, from every day experience we know 
that we can generate a heat not only by transferring en- 
ergy from say hot to cold object but by doing some me- 
chanical work. For example, if we take a container with 
a gas isolated from the environment and then perform 
some cyclic process over this gas jT(|, e.g. compress and 
decompress it, then we will generate the heat unless we 
do it very slowly (adiabatically) . In quantum mechanics 
for a cyclic process the heat will be equal to the energy in- 
crease caused by the transitions between different energy 
levels and thus it can be microscopically defined even for 
an isolated system. If two systems are brought to a con- 
tact for certain time and then again separated then heat 
can be defined in a similar fashion: energy change of a 
system due to transitions between different energy lev- 
els. On the other hand by work in thermodynamics one 
usually understands reversible part of the energy change, 
which is due to the change in the energy level structure 
but not due to the transitions. In a closed system then 
the total energy at each moment of time is 

£(t)=J2 E n(t)pnn(t), (1) 
n 

where p nn (i) are the diagonal elements of the density ma- 
trix in the instantaneous adiabatic basis. This expression 
can be rewritten as 

£{t) =J2 E n(t)Pnn(Q) + ^ E n (*) (pnn (*) " ft»„(0)) 
n n 

= e a (t) + Q(t), (2) 

where p nn (0) are the diagonal matrix elements of the ini- 
tial density matrix. The first term corresponds to the 
energy in the adiabatic process. Indeed as we know from 
quantum mechanics in the adiabatic limit there are no 
transitions between energy levels and thus diagonal ele- 
ments of the density matrix do not change. The second 
term in Eq. ((2]) is related to the transitions between dif- 
ferent levels and thus corresponds to the heat. Of course 
in large systems it is impossible to completely avoid tran- 



sitions between energy levels However, for slow pro- 
cesses the amount of states available for such transitions 
is small. Thus their contribution to the energy is negli- 
gible and we recover the adiabatic limit. There is an- 
other reason why splitting of energy in Eq. ([2]) into the 
adiabatic part and the heat is consistent with thermo- 
dynamics. The first (adiabatic) term in this expression 
depends only on the state of the system, i.e. only on the 
initial state of the system and on the instantaneous pa- 
rameters of the Hamiltonian. At the same time the heat 
is clearly a function of the process because it depends on 
the transitions between different energy levels during the 
evolution. 

The heat, in turn, is intimately connected to the second 
law of thermodynamics (or simply the second law) , which 
is the fundamental law of physics. In the most com- 
mon formulation it states that the entropy of an isolated 
system can either increase or stay constant. However, 
there are other thermodynamically equivalent formula- 
tions of the second law [3(. In particular, the Kelvin's 
(Thompson's) formulation, used e.g. to prove impossi- 
bility of the perpetuum mobile of the second kind, states 
that it is impossible to create work by extracting heat 
from an isolated system. This statement implies that if 
we perform a cyclic process over some system then we 
can only increase its energy. Thus the heat in Eq. @ 
should be nonnegative if we start from thermal equilib- 
rium. And indeed this statement was microscopically 
proven m Refs. HQ. In Ref. Q it was further shown 
that Q(t) > for non-cyclic processes (and still initial 
thermal equilibrium state) as long as there are no un- 
avoided level crossings between different energy levels. 
This statement is called the minimum work principle. In 
Ref. 0] a more general statement was formulated that 
the work on the system or equivalently heat is nonnega- 
tive for any cyclic process if the initial density matrix is 
passive, i.e. i) it is diagonal and ii) microscopic occupan- 
cies of different energy levels are nonincreasing functions 
of energy, i.e. (p n - p m ){E m - E n ) > 0. 

The main purpose of this paper is to find a microscopic 
expression for the heat in the system under the assump- 
tions and relate it to many-body transition probabilities. 
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The other purpose of this work is to define the frame- 
work for the perturbation theory around the adiabatic 
limit. For this purpose it is natural to work in the adi- 
abatic or instantaneous basis. This basis diaogonalizes 
the Hamiltonian at each moment of time and thus if the 
Hamiltonian changes in time so does the basis. In Ref. Q 
it was already demonstrated how one can use this basis 
to calculate scaling properties of the density of excita- 
tions near critical point doing the first order perturbation 
theory in the rate of change of external parameters. In 
Ref. [9( it was shown how using this perturbation theory 
in the adiabatic basis for a particular sine-Gordon model 
one can find the number of excitations and the heating 
using form-factors. Basically finding dynamical observ- 
ables was reduced in that work to taking integrals over 
known static matrix elements. 

Here we find a microscopic expression for the heat and 
more generally for diagonal matrix elements of the den- 
sity matrix under the conditions that the initial density 
matrix is stationary (diagonal): 

Pnn(t)=P° n + J2(P°m-Pn)\Inm(t)\ 2 , (3) 

m 

where p° is the initial value of the occupancy of the n-th 
energy level, I n m{t) = Inm are the microscopically de- 
fined transition amplitudes, which will be explicitly de- 
fined below. Remarkably the sum in Eq. ([3]) contains 
only squares of the absolute values of I nrn suggesting 
that | Inm | 2 can be interpreted as a transition probabil- 
ity from m-th to n-th state. Note that since \I nn \ 2 in 
Eq. ([3]) is arbitrary one can always choose it such that 
X^m l-^nml 2 = 1< As we point below the unitarity of the 
evolution requires that transition probabilities satisfy the 
sum rule: 

^ ^ |-^mn| — ^ ^ |-^nm| ■ (^) 

rn rn 

In general \I mn \ ^ \I nm \ though the transition probabil- 
ities are equal to each other within (i) the Fermi-Golden 
rule, which applies, for example, for infinitesimally short 
transition times, or (ii) for cyclic processes where the ex- 
ternal parameter symmetrically changes in time. Using 
Eq. ([3|) to find the total energy of the system and sub- 
tracting the adiabatic part we find 

Q = J! En {Pm - Pn)\In7n\ 2 = Pm( E n ~ E m )\I nm \ 2 . 

m.n m,n 

(5) 

For a cyclic process where final and initial energies are 
identical: E n = E® one can show Q(t) > if the initial 
state is passive [7| . If the transition probabilities satisfy 
detailed balance |/„ m | 2 = |I nm | 2 then the proof trivially 
follows from the symmetrization of the sum with respect 
to the indices m and n. In general the proof is more 
subtle but the statement is still correct [1, @] ■ From the 



explicit expressions below we find that in the adiabatic 
limit all I nm — > as long as E n ^ E m and we recover 
Q — > 0. For non-cyclic process the heat is still positive 
under the same assumptions as long as there are no un- 
avoided level crossings [6|. 

Let us now make a few remarks on Eqs. ([3]) and ([5]): 

• In derivation of Eq. ([3]) it is important that the 
initial density matrix is diagonal only at the ini- 
tial time, the following evolution can be arbitrary 
and the density matrix at any later time t is not 
necessarily diagonal. 

• The Gibbs distribution p n oc exp(— pE n ) is obvi- 
ously a passive state resulting in non- negative heat. 
However, the same statement equally applies to ini- 
tial generalized Gibbs distribution functions sug- 
gested for integrable systems [l^ . 

• The corollary of the Kelvin's statement is also 
true: if one starts from an overheated configuration 
where higher energy states have higher occupancy 
then the heat after a cyclic will be non positive, i.e. 
it will not be possible to do a positive work on the 
overheated system. 

• The state where all energy levels are equally popu- 
lated (and the initial density matrix is stationary) 
is neutral in a sense that its diagonal matrix ele- 
ments can not be changed in any process. Corre- 
spondingly one can not add or extract the energy 
from this state. 

In the remainder of the paper we will give the details 
of the derivation of Eq. ([3]). First we will show how the 
derivation works for a two-level system and then we will 
generalize the result. 

Two-level system. Suppose we are dealing with a 
Hamiltonian two level system. We will work with the 
Schrodinger equation for the wave function describing the 
system and then take the average over the initial density 
matrix. Alternatively one can directly work with the Li- 
ouville equation for the density matrix. In the adiabatic 
basis (co-moving with the Hamiltonian) the Schrdinger 
equation for the wave function takes the form [8|: 

-± = -R2,i(t)ai(t), = -Ri, 2 (t)a 2 (t), (6) 

where a\{i) and a 2 (t) are the coefficients of the expan- 
sion of the wave function in the adiabatic basis describ- 
ing our two level system ip(t) = ai(t)\l) + a%{t)\2). The 
states |1) and |2) are the eigenstates of the instantaneous 
Hamiltonian: H{t)\n) = E n (t)\n), n = 1,2. Finally 
R2,i(t) = -iZ* j2 (t) is given by 

R 2 ,i{t) = (2|9 t |l)exp (ij (E 2 (f) - E^dt'^j . (7) 
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For simplicity we omitted the Berry phase term here. 
The latter can be always reabsorbed into definition of 
i?2,i and thus it is not going to affect our conclusions. 
We work in units where H = 1. The condition i? 2 ,i = 
— R* 2 ensures the unitarity of the evolution of the wave 
function. Eqs. © can be rewritten in the integral form: 

aa(*)=oa(0)- f R 2 . 1 (t / )a 1 {t')dt' (8) 
Jo 

oi(t) = Oi(0) + [ (t')dt' (9) 

Jo 

We are interested in the diagonal elements of the den- 
sity matrix at time t (by analogy one can also-derive 
off-diagonal matrix elements). For example we have 
P22.it) = 0-2 (i) | 2 , where the overline implies averaging 
over the initial density matrix. A similar expression can 
be found for pu(t). Using Eqs. §8§ and §§§ we find 



M*)=P°+ / dt'dt"R 2A (t')Rl 1 (t' , )a 1 (t>)a*(t>' 



dt' a*{Q)R 2 ,i{t')ai{t') + a 2 (0)R* 1 (t , )a* 1 (t l ). (10) 



Let us first prove our statement assuming that we are 
close to the adiabatic limit and the system is only slightly 



excited [17j. Using the fact that 02(0)01 (0) = we find 

P22 (t) ~p° 2 -p° 2 [G 2 .i o G* 2A + (?2,i o G 2) i] 
+ P° 1 \G 2 .i\ 2 = P 2 + (Pi-P 2 )\G 2 ,i\ 2 . (11) 

Here G„. m and their convolutions are defined as follows: 

G 2 ,i = G 2A (t) = f i? 2 ,iW, 
Jo 

G2,i°G^= f dt'R 2A {t') f dt"R 2A (t"). (12) 
Jo Jo 



Note that these convolutions can be also expressed 
through time-ordered products. To get from the first 
to the second line in Eq. (fTT|) we used an obvious iden- 
tity G 2i i o G 2 1 + G 2l o G 2i i = |G 2 ,i| 2 . In the adiabatic 
limit (assuming that the system is not degenerate) we 
have G2,i = and hence we have p 22 {t) — p\. In general 
G 2 .i ^ if the transitions between the two levels are not 
forbidden and as a result p 22 (t) > p 2 as long as Pi > p 2 , 
which is the condition of passivity for a two-level system. 

Beyond the perturbation theory we can solve Eqs. © 
and ([5]) by iterations and substitute the solution into 
Eq. pop . Since we give a general derivation for a multi- 
level system below we will only state the final result: 



P2 2 (t)= P ° 2 + (Pl-P 2 )\l2,l\ 

T G 2 .l 
^2,1 = 



(13) 

The quantity I 2 ,i is understood in terms of the Taylor 
expansion: 7 2 ,i = G24 — G 2 ,i ° G%^ o G 2 .\ + . . . Alterna- 
tively I 2 ,i can be represented as a time ordered exponent 
(see Eq. ([20]) below. 

Multi level generalization. We can repeat now a sim- 
ilar derivation for a general system containing multiple 
energy levels. We again will work in the adiabatic basis. 
Then instead of Eq. ((8|) we find 



a n {t) = a„(0) - 



R n , m {t')a m (t')dt' , (15) 



where the prime over the summation index implies that 
the term with m = n is excluded. For simplicity we 
again assume that there is no Berry phase appearing in 
the problem, the latter can be always reabsorbed into the 
definition of R n ,m- Using Eq. (|15[) one immediately finds 
the generalization of Eq. (HI 



pnn{t) =p° + 



m >P n 



dt'dt"R n , m {t')R* n>p (t")a m (t')a; (t") 



dt' al{0)R n ,m{t')a m {t') + a n {0)R^ m {t')a* m {t>). 



(16) 



Let us again assume that we are close to the adia- 
batic limit and do the perturbative expansion in G n _ m — 

Jo* Rn.m{t')dt' . In the leading order we find 

Pnn(t)=p°+J2(P°m~Pn)\G m ,n\ 2 . (17) 

m 

In the adiabatic limit G n ^ m — > for all n ^ m as long as 



E n ^ E m and we recover p nn (t) = Pn- If the process is 
not adiabatic and the transitions between different levels 
are not forbidden then diagonal matrix elements of the 
density matrix change in time. The equation above can 
be generalized beyond the first order of the perturbation 
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expansion. We can formally solve Eq. (1 1 5j) : 
a n (t) = <x n (0) - ^ ln, m am(0), 



where 



G 



l + G 



(18) 



(19) 



Note that formally I nm can be expressed through a time 
ordered exponent 



I = — T T exp 



Jo 



(20) 



where 1 is the unity operator. Next we substitute this 
expression to Eq. p^|) to get: 



(21) 



Now instead of checking the desired properties of the 
Green's functions leading to Eq. ([3]) we will use the fol- 
lowing trick. Consider the normalization of the wave 
function ^2 n \a n (t)\ 2 = 1. This implies 



1 = 1 + E Pn l J ™ I' " E Pn + 7 n, 



(22) 



Clearly the normalization should be conserved for any 
choice of p n , therefore 



Inn + I n n ~ ^ t \Imn\ — ^ ] |Jri 



(23) 



The identity above actually follows from the unitarity of 
the evolution and it is equivalent to the optical theorem 
in scattering theory [l5|. Note that in the sum above 
we are allowed to interchange indices m and n because 
(I'I)nn — {H^)nn- This identity trivially follows from 
the fact that / is the difference between unity and uni- 
tary operators (see Eq. ([20]) ). Substituting Eq. (|23|) into 
Eq. (J2l]) we recover Eq. (J3|). 

The result ([3|) can be used for finding expectation val- 
ues of other observables commuting with the Hamiltonian 
or of the time averaged values of arbitrary observables 
following a quench. If is such an observable then 

(H) = J2 VnnPn + E "™ Oh ~ PnWnm (t)\ 2 , (24) 
n rn,n 

where overline implies the time averaging, which is nec- 
essary if does not commute with the Hamiltonian. 
As long as the matrix elements of ft satisfy the condi- 
tions of passivity: for p n > p m we have tt nn > Q mm 
(0„„ < Slmm) the the expectation value of this opera- 
tor can either decrease (increase) in an arbitrary cyclic 
process. In particular, one can choose A to be equal 



to any power of the Hamiltonian an make similar state- 
ments about the moments of the Hamiltonian. We also 
note that Eq. (|24| gives expectation values of steady state 
values of arbitrary stationary observables not commuting 
with the Hamiltonian. This is true because if the Hamil- 
tonian does not change in time then all off-diagonal el- 
ements of the density matrix oscillate in time and thus 
average to zero once the system reaches the steady state. 

The expansion (fT9|) together with Eq. ([3]) can be used 
to develop perturbation theory in the quench rate 5 (or 
the rate of change of external parameters) near the adi- 
abatic limit. In Refs. [2j, |8j, [9j such expansion in the 
leading order was used to analyze slow dynamics for par- 
ticular situations. It was shown that this perturbation 
theory can be applied when the conventional perturba- 
tion theory in interaction strength diverges, for example, 
in low dimensional systems particularly near phase tran- 
sitions. We also emphasize that Eq. ([3]) was obtained 
under the assumption of the initially diagonal density 
matrix. If this is not true then the expression for p nn (t) 
will acquire additional terms linear in Pn,7n 

(0) for n 7^ m. 

Those terms usually strongly fluctuate in the phase space 
in complex systems. Then these additional terms will av- 
erage to zero and thus will not affect the Eq. ([3]). 

In conclusion we derived a microscopic expression for 
the diagonal elements of the density matrix ([3]) in the 
instantaneous adiabatic basis provided that initial den- 
sity matrix is stationary (diagonal). This expression can 
be used to find explicit formulas for the heat, which 
is consistent with the Kelvin's formulation of the sec- 
ond law of thermodynamics. In the adiabatic basis all 
transition probabilities are expressed as squares of the 
matrix elements of the many-body evolution operator. 
The time only enters as an external parameter so there 
is no qualitative difference between time dependent and 
time-independent problems (if we are interested only in 
diagonal matrix elements of the density matrix). The 
Eq. ([3]) can be used to construct an adiabatic pertur- 
bation theory, to perform Monte-Carlo simulations for 
time-dependent problems, and to derive general connec- 
tions between various stationary observables in the out 
of equilibrium systems. 
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